U{1) Slave-spin theory and its application to Mott transition in a multi-orbital model 

for iron pnictides 
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A 17(1) slave-spin representation is introduced for multi-orbital Hubbard models. As with the Z2 
form of L. de'Medici et al. (Ref. [J), this approach represents a physical electron operator as the 
product of a slave spin and an auxiliary fermion operator. For non-degenerate multi-orbital models, 
our f/(l) approach is advantageous in that it captures the non-interacting limit at the mean-field 
level. For systems with either a single orbital or degenerate multiple orbitals, the U{1) and Z2 
slave-spin approachs yield the same results in the slave-spin-condensed phase. In general, the U{1) 
slave-spin approach contains a U{1) gauge redundancy, and properly describes a Mott insulating 
phase. We apply the U{1) slave-spin approach to study the metal-to-insulator transition in a five- 
orbital model for parent iron pnictides. We demonstrate a Mott transition as a function of the 
interactions in this model. The nature of the Mott insulating state is influenced by the interplay 
between the Hund's rule coupling and crystal field splittings. In the metallic phase, when the Hund's 
rule coupling is beyond a threshold, there is a crossover from a weakly correlated metal to a strongly 
correlated one, through which the quasiparticle speactral weight rapidly drops. The existence of 
such a strongly correlated metallic phase supports the incipient Mott picture of the parent iron 
pnictides. In the parameter regime for this phase and in the vicinity of the Mott transition, we find 
that an orbital selective Mott state has nearly as competitive a ground state energy. 



I. INTRODUCTION 

Many important questions remain on the physics of 
the iron pnictides and related iron-based high- Tc super- 
conductors^ri. One central issue is the strength of elec- 
tron correlations these systems contain. The metallic na- 
ture and the coUinear antiferromagnetic (AFM) ground 
stated in the parent pnictides can arise within a weak- 
coupling approach, in which the Fermi surface nesting 
plays an important role^. On the other hand, the large 
room-temperature electrical resistivity, showing the bad 
metal behavior with the electron mean free path on the 
order of the inter-electron spacing, as well as the suppres- 
sion of Drude weight in optical conductivityi2r22, the well 
defined zone-boundary spin wavesi^, and the renormal- 
ization of LDA bandstructure in the ARPES mcasure- 
mentsii provide evidence for sufficiently strong electron 
correlations such that the system is in close proximity 
to a Mott transition with dominant contributions to the 
spin spectral weight from quasi-local momcntai^^— . 

This incipient Mott picture is further supported by the 
properties of iron chalcogenides^"— . In either the 11- or 
the 122-chalcogenides, the magnetic ordering wave vector 
and the large magnetic moment^"— (ranging from 2 to 
3.4 ^b) can hardly be explained by a Fermi surface nest- 
ing mechanism, but is readily understood within a quasi- 
local moment model^"— . It has also been shown that the 
band narrowing effect, either due to the expansion of the 
iron lattice unit cell in a iron oxychalcogenide^, or from 
the ordered iron vacancies in 122 iron selenides^i^i^, 
may drive the system through the Mott transition to a 
Mott insulator. 

In general, the degree of electron correlations can be 
measured by the ratio U/D, where U refers to a charac- 
teristic interaction strength, such as the Coulomb repul- 
sion in a Hubbard model, and D is the full bandwidth, a 



scale of the kinetic energy of the system. A Mott tran- 
sition separating the metallic and the Mott insulating 
phases takes place at Uc ~ D. In the metallic phase, 
electron correlations can alternatively be measured by 
the quasiparticle spectral weight Z, which is unity at the 
non-interacting limit U = 0, and vanishes at the Mott 
transition U = Uc- The incipient Mott picture relies on 
the existence of a (putative) Mott transition, and assume 
that the system is not too far from this transition so that 
Z is relatively small. 

The above considerations suggest that it is very im- 
portant to theoretically investigate how a metal to Mott 
insulator transition takes place in a model which is appli- 
cable to the parent iron pnictides and/or chalcogenides. 
In a previous paper, by studying a two-orbital and a 
four-orbital model, the authors have shown that a Mott 
transition generally exists in these models^. But both 
models are at half-filling. By contrast, in the parent iron 
pnictides, six electrons occupy five 3d orbitals of each Fe. 
This means that the system is away from half-filling. It 
would be important to investigate extent to which the 
Mott insulating states persist in such a situation. 

Historically, the Mott transition in a single-orbital 
Hubbard model has been studied by using various tech- 
niques, such as Dynamical mean-field theory (DMFT) 
and Gutzwiller approximation. Among these methods, 
the Kotliar-Ruckenstein slave-boson method^ has been 
broadly used. In this theory, a Mott transition at finite 
U /D is obtained already at the mean-field level. How- 
ever, it is difficult to apply this approach to multi-orbital 
systems as it would introduce 4^^ slave boson fields for 
a model with M orbitals; the number of variational vari- 
ables is already huge even for M = 5. Recently several 
other slave particle theories following the idea of charge- 
spin separation have been proposed. These include the 
slave-rotor theory^, and the slave-spin theory^. In both 
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theories, a slave bosonic variable (a quantum 0(2) rotor 
in the slave-rotor theory, and a quantum 5=1/2 spin in 
the slave-spin theory) is introduced to carry the electric 
charge, and an auxiliary fermion (the spinon) carries the 
spin of the electron. The metallic phase corresponds to 
the state that the slave particles are Bose condensed, so 
that charge excitations arc gaplcss along with the spin 
excitations. By contrast, the Mott insulator corresponds 
to the state that the slave particles are disordered and 
gapped; charge excitations are gapped while the spin ex- 
citations remain gapless. The slave rotor method is very 
economical because it introduces only one rotor per site 
taking account of the total charge. It is also very efficient 
if the interactions have an SU{2M) symmetry for an M- 
orbital system. However, it can not be easily applied 
to systems with a non-zero Hund's rule coupling which 
breaks the SU{2M) symmetry. It is also not convenient 
to handle systems that exhibit strong orbital dependence, 
such as the orbital selective Mott transition (OSMT). 
The slave-spin theory overcomes these drawbacks by in- 
troducing a slave spin for each orbital and spin flavor. 
Compared with the Kotliar-Ruckenstein slave-boson the- 
ory, it is still very economical because it introduces only 
2M slave spins per site. It has been successfully used 
to study the Mott transition and OSMT in multi-orbital 
systems with a non-zero Hund's rule coupling^. There 
are however several issues with the slave-spin represen- 
tation. First, it has difficulties when applied to a multi- 
orbital system in which one or more orbitals are away 
from half-filling due to crystal field splitting. Secondly, 
in its original construction, the slave-spin representation 
has a Z2 gauge redundancy. This makes the spinous carry 
both spin and charge currents, and causes difficulties in 
describing a Mott insulating phase — . 

In this paper we propose an XY slave-spin theory that 
is free of these issues. This slave-spin theory has a C/(l) 
gauge redundancy, and properly describes Mott insulat- 
ing phases. We develop a mean-field theory, which can be 
applied to multi-orbital systems with a non-zero crystal 
field splitting. For a model with a single orbital or degn- 
erate multiorbitals, on the other hand, our U{1) slave- 
spin mean-field theory and the Z2 slave-spin theory^ 
give the same results in their spin-condensed phases. We 
then apply our formulation to study the Mott transi- 
tion in a five-orbital model for the parent iron pnictides. 
We establish the existence of a Mott transition in this 
model. Both the nature of the metallic and Mott in- 
sulating phases are strongly affected by the interplay of 
Hund's rule coupling and the crystal field splitting. A 
crossover to a strongly correlated metallic state exists 
when the Hund's rule coupling is beyond a threshold. 
The existence of this state is in agreement with the in- 
cipient Mott picture. 

The rest of the paper is organized as follows. In Sec.HIl 
we first introduce our construction of the U{1) slave-spin 
theory, and develop a mean-field theory based on this new 
construction. We also compare our construction with 
the slave-rotor and Z2 slave-spin theories. In Sec. IIIII we 



apply the (7(1) slave-spin mean-field theory to study the 
Mott transition in a five-orbital model for the parent iron 
pnictides, and show how the transition is affected by the 
interplay between Hund's rule coupling and the crystal 
field splitting. Finally Sec. IIVI contains some concluding 
remarks. 



II. METHOD 

In this paper we arc interested in the metal-to-insulator 
transition in a multi-orbital system, for which the Hamil- 
tonian reads 



H — Hf) + Hint ■ 



(1) 



Hq contains the tight-binding parameters among the 
multiple orbitals, 

^0 = ^ E tfdl^d.p, + E(^" - Mcad^ca, (2) 
ijaj3(T iaa 

where dj^^ creates an electron in orbital a with spin a at 
site i, Aq, is the on-site energy reflecting the crystal held 
splitting, and n is the chemical potential. Hint contains 
on-site Hubbard interactions 



Hi, 



U 



E 



i.a</3.<T 



-J{d\aAasd\p^d^pa " d\^j\^^d,fi„d,fis)^ (3) 

where riiacr = dj^^diaa- In this model, U, U', and J re- 
spectively denote the intraorbital repulsion, the interor- 
bital repulsion, and the Hund's rule exchange coulping. 
In the following, we wiU take U' = U — 2 J4i 



A. U{1) slave-spin theory 

To study this multi-orbital model, we implement the 
idea of charge and spin separation of the d electrons by 
using the slave-spin approach. For each orbital and spin 
flavor we rewrite the electron creation operator to be 



ick.a iaa J iaa "> 



(4) 



where S^^^ is a ladder operator of the slave quantum 
5 = 1/2 spin carrying the charge degree of freedom of 
the electron, and fj^^ is a spinon creation operator. We 
further enforce a constraint for each site 

^ia<T ~ f iaa- f iaa ^ — , (5) 

which restricts the Hilbert space to the physical one. 
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Note that our slave-spin formulation in Eg . (|4]) is dif- 
ferent from that introduced in Refs. [l] andl48l in which 



where 



(6) 



(7) 



with daa being a complex number. In that formula- 
tion, the gauge redundancy is reduced from U{1) to Z2 
in Eq. (|6]) and Eq. ^ due to the mixing between S^^^ 
and 5*",^. (hence they are referred to as a Z2 slave-spin 
representation) ; as a consequence, the slave spins can not 
carry the U{1) charge. 

In our Eqs. (H]),©, there is a U{1) gauge redun- 
dancy corresponding to fLae'''^"'"' and 5'+^ ^ 
S'j^g.e'^'"" . The slave spins carry the C/(l) charge, sim- 
ilarly as the slave rotors — . We will refer to this as a 
U{1) slave-spin theory. In our construction, the phase 
that corresponds to disordered slave spins (preserving the 
U{1) symmetry) but with gapless spinous corresponds to 
a Mott insulator. 

Next we develop a mean-field theory based on the con- 
struction of Eq. (|4]) and Eq. ([5]) . A naive mean- field the- 
ory based on Eq. ^ would not produce the correct quasi- 
particle spectral weight in the non-interacting limit^. 
To make progress, we rewrite the slave spin operators in 



their Schwinger boson representation: S^^^ 



s;. 



and Sf, 



The constraint in Eq. ([S]) then becomes 



h'^ h- 



^fiacrfiacr — 1- Here wc need to introduce 
an extra constraint al^^Uiaa + bl^^biaa = 1 so that the 
Schwinger bosons represent S = 1/2 spins. We then see 
that flictcr and biaa are hard-core bosons. In light of the 
Kotliar-Ruckenstein slave-boson mean-field theory^, we 
now define a dressed operator in the Schwinger boson rep- 
resentation which automatically takes account this con- 
straint: 



P+ b 



(8) 



where P,, 



1/^1/2 + 6 ±{al 



bLab^o..)/2, 



± 



and S is an infinitesimal positive number to regulate 

and a^b are equivalent in the physical Hilbert space. 
In the Schwinger boson representation, Eq. becomes 



iota iaa J tacr ' 



(9) 



At the mean-field level, we treat the constraint Eq. ([5]) 
on average by introducing a Lagrange multiplier, and de- 
compose the boson and spinon operators. We obtain two 
mean-field Hamiltonians respectively for the spinous and 



the Schwinger bosons: 



— ^ ^ ^if {^iaa^3l3cr)fiaafjl3<T 



+ ^(A„ - A,,, - ^l)flJ^c.a, (10) 

icxa 

rrmf \ ^ ,cxj3 i p-\ p \ t 

+ T.^(f^^c.a-nl^)+HL, (11) 



where 



• ) denotes the mean-field value, nf^^ = 
'^iaa'^iatTi and Xiaa IS the Lagrange multiplier to han- 
dle the constraint in Eq. ([5]). The quasiparticle spectral 
weight is defined as Ziaa = \{ziaa)\'^- In Eq. (fTT|) H^^ 
refers to the interaction Hamiltonian in the Schwinger 
boson representation. It can be obtained by rewriting 
Eq. ([3]) in the slave-spin representation i/jnt -ffint(S),— 
then substitute the Schwinger bosons for the spin oper- 
ators. The mean-field Hamiltonian i?™^ has an internal 
U{1) symmetry of the bosons. For a single orbital, it is 
a Bose Hubbard model for two species of bosons, and is 
equivalent to a model of interacting XY spins in a mag- 
netic field. At commensurate fillings, by breaking the 
internal U{1) symmetry, this model has a phase transi- 
tion from a bosonic Mott insulator to a superfiuid with 
decreasing the interactions. These two phases correspond 
to the Mott insulating and metallic states in the original d 
electron problem. We then approach the Mott transition 
from the ordered phase, and further adopt a single-site 
approximation to Eq. ([TU]) and Eq. ([TT|) with the decou- 
pling 4aa^i/3<T ~ {4acr)Zjf3a + zja^iZjfSa) - {zLa) i^J 0<t) ■ 

For simplicity, we focus on the paramagnetic phase, and 
assume translational symmetry. These allow us to drop 
the spin and site indices in the formulas. The mean-field 
boson Hamiltonian then reads 

a/3 



(12) 



where e'^^ = E^.jtfifLfjP^)/^- In Eq.^dHD, we 
Taylor-expand Zq, and z^ in terms oi A — (A) (where 
A — fi"' , n}' , a'' b) , and keep up to the linear terms in 
A — (A) , obtaining 

4 « 4 + {sDvaiK - < - (2^4 - 1)], (13) 

where 4 = (P+)at &„(p-), = (2n/-l)/[4n/(l-n/)], 

'4 = ^ T,kifLfkc.)^ and ni = {n'^) = 1- (n^) from the 
constraints. We find that Eq. already gives good 
mean-field results, and will hence drop the higher order 
terms in the expansion. Note that the approximate form 
of in Eq. is only used to simplify H™^ , but can 
not be fed into Eq. © to calculate the electron Green 
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functions since the approximate operator behaves differ- 
ently from the original one in the physical Hilbert space. 
Nevertheless, = {za). With Eq. Eq. ^ is 

then approximated to be 



m 



afi 



E 



+ e^Vo. ] (K - fi'a) + Htnt, (14) 



where Cq, = J^pi^"'^ {^a){^l3) + c.c). Further using the 
constraint Eq. ([5]) , we can move the term proportional to 
r]a to HJ^^ by introducing an effective on-site potential 
fla = 26"?]^. The resulting mean- field Hamiltonians are 
then 



the spinon bandstructure is identical to the tight-binding 
one, and the non-interacting limit is properly recovered. 

However, we find that at the mean-field level and in 
the symmetry broken phases of the bosons/spins, the two 
theories have very similar forms. To see this explicitly, we 
compare the mean-field Hamiltonians of the two theories 
in the slave-spin representation. In this representation, 
the mean- field Hamiltonian of the U{1) slave-spin theory 
can be obtained by performing a Schwinger-boson-to-spin 
mapping to the Hamiltonian Hf^ in Eq. (dH). Hf^ then 
reads 



E 

ai3 



Hf 



where 



where e^*^ = wJ^ij ^if ^-nd Sap is Kronecker's 

delta function. Eq. ([T5|) and Eq. (fT6|) are the main for- 
mulation of our slave-spin mean-field theory. The mean- 
field parameters (zq,) and Aq, can then be solved sclf- 
consistently. In the non- interacting limit, it is easy to 
check that = = 1 can be achieved by taking 

Act = Act; the quasiparticle weights are equal to 1 as it 
should be. 



B. Comparison with the Z2 slave-spin theory 

Here we compare our U{1) slave-spin theory with the 
Z2 slave-spin theory. One advantage of the U{1) slave- 
spin theory over the Z2 slave-spin approach is that it can 
be directly generalized to the multi-orbital systems with 
non-zero crystal field splitting and/or away from half- 
filling. As an example, wc consider a two-orbital model 
at half-filling and with equal bandwidth but with a finite 
crystal field splitting. In the non-interacting limit, we 
expect the quasiparticle spectral weight ^ct=i.2 = 1 and 
the spinon bandstructure is identical to the tight-binding 
dispersion of the d electrons. The Z2 slave-spin mean- 
field theory failed to obtain these results. This is easy to 
see: Both of the two orbitals are away from half-filling 
in the presence of the crystal field splitting. According 
to Eq. p7|) . to obtain Zi(^2) ~ |(0i(2))P = 1, it is nec- 
essary that Ai = — A2 ^ 0. In absence of the potential 
p, in Eq. ([15]), this already distorts the bandstructure of 
the spinous from the original tight-binding form. Hence 
one can not obtain the desired spinon filling as required 
by the constraint. However, in our theory, Za = 1 is 
guaranteed by the condition Aq = These two poten- 
tials cancel out exactly as seen in Eq. ([TSl) . Therefore, 



+i?i„t(S), 



{Op)Oi)+5apKSl] 



(15) 



(17) 



(18) 



(16) 



P± = ^1/2 + 5 ± Si, and = |(0,)|2. Surprisingly, 
we see that the mean-field Hamiltonian of the Z2 theory 
takes exactly the same form as in Eq. (|17p if define 



01 



(19) 



Interestingly, the two definitions in Eq. (fT8|) and Eq. (fT9)) 
gives the same quasiparticle spectral weight Za ~ 
I (Oct) p. In models with a single orbital or degenerate 
multiple orbitals, (la in Eq. ([T^ becomes orbital inde- 
pendent and can thus be absorbed into the chemical po- 
tential. Therefore, in these cases, in the metallic phase, 
the U{1) slave-spin mean-field theory and the Z2 theory 
gives the same results (up to a constant in the free ener- 
gies). It should be stressed that, even for single-orbital or 
degenerate multi-orbital models, the agreement between 
the two theories is obtained only in the ordered phase of 
the slave spins. 

Generally, the constructions in the two formulations 
are different in the sense already mentioned. In the U{1) 
slave-spin theory, the operator equation Eq. (|4]) has a 
C/(l) gauge redundancy, and allows a proper Mott insu- 
lator phase. 



III. MOTT TRANSITION IN A FIVE-ORBITAL 
MODEL FOR IRON PNICTIDES 

In a previous paper— 1^ we have discussed the metal- 
to-insulator transition in two- and four-orbital models. 
For both models wc find a metal-to-Mott-insulator tran- 
sition at finite U. We also find that the nature of the 
phase transition and the critical value Uc can be strongly 
influenced by the Hund's rule coupling. These results are 
based on the assumption that the system under study is 
at half-filling. But the parent iron pnictides are away 
from half-filling, with six d electrons occupying five or- 
bitals in each iron atom. To consider this more involved 
case, here wc use the U{1) slave-spin mean-field theory to 
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FIG. 1. Ground-state phase diagram of the five-orbital model 
at n = 6.0. The solid curve with symbols shows the boundary 
between the metallic and Mott insulating phases. The dotted 
line shows a crossover in the metallic phase where the system 
changes from a weakly correlated metal to a strongly corre- 
lated metal. The dashed line indicates a low-spin to high-spin 
transition in the Mott insulating phase. The atomic config- 
urations corresponding to the low-spin and high-spin Mott 
states are illustrated on the right side. 



study the five-orbital Hubbard model for iron pnicitdes. 
We take the tight-binding parameters as those proposed 
in Ref. for the parent LaOFeAs. The interaction part 
of the Hamiltonian is as given in Eq. Q. For simplic- 
ity, here we consider only the density-density interactions 
and drop the spin-flip and pair- hopping terms in Eq. . 
The results with the full interactions are qualitatively 
similar—. 

Fig. [T] shows the ground-state phase diagram at elec- 
tron filling n = 6 of the five-orbital model in the J- 
U plane. We find a single transition from a metal to 
an insulator. At J = 0, J7c « 13.1 eV. In this case, 
we have also independently determined U^^ w 11.0 eV 
from the slave-rotor mean-field theory^. For a model 
with M-fold-degenerate orbitals, the ratio Uc/U^^ = 
{M + 1)/MJu^ This relation approximately holds in the 
five-orbital model, in which we find Uc/U^^ w 1.19. It 
can be understood by the fact that in this model, the 
largest crystal field splitting (A w 0.5 eV) is relatively 
small compared to the full bandwidth {D « 4.0 eV), so 
the orbitals are nearly degenerate. 

A non-zero Hund's rule coupling strongly affects the 
MIT and the nature of both the insulating and metallic 
states. In the insulating phase, we find that the degen- 
erate xz and yz orbitals are always at half filling, and 
hence in a Mott insulating state; while the — or- 
bital is fully occupied, in a band insulating state. The 
transition is therefore an orbital selective MIT.— Due to 
the interplay of J and crystal field splitting A, the other 
two orbitals can be either in a Mott insulating ( J > A), 
or in a band insulating state (J < A). Accordingly, the 




U(eV) 



FIG. 2. (a): Evolution of the quasiparticle spectral weight 
with U a.t n — 6.0 and J/U = 0.1. (b): Same as (a) but at 
J/U — 0.25. The inset shows the difference of free energies 
(A/) between three competing states for the same model pa- 
rameters, with /m, /mi, and fos respectively denote the free 
energies of the metallic, Mott insulating, and orbital selective 
Mott states. 



five-orbital model can accommodate either a high-spin 
{S = 2) or a low-spin (S = 1) Mott state, as illustrated 
in Fig. [1] These two states are separated by a low-spin- 
to-higli-spin transition inside the insulating phase. 

In the models at half filling, Uc decreases monotoni- 
cally with increasing J/U. But the phase diagram of the 
five-orbital model shows a significant difference: as J/U 
is increased from zero, Uc first decreases for J/U < 0.1, 
but then increases with J/U for J/U > 0.1. We can 
understand this by estimating and comparing the Mott 
gaps in the four- and five-orbital models in two limiting 
cases: < J <C A and J > A. The Mott gap G'^ mea- 
sures the distance between the upper and lower Hubbard 
bands in a Mott insulator. It can be approximated by 
G^^ G^ ~ D; here G"^ = + E''-^ - 2£:" is the 

excitation gap in the atomic limit, and E"' is the energy 
of the atomic configuration with n electrons. A reason- 
able estimate of Uc can be obtained from G^ ^ D. For 
< J ^ A, the configuration for the undoped com- 
pound is the S — \ low-spin state. For this configu- 
ration, G^ — U + J UY both the four- and five-orbital 
models. Hence when J/U is small, in both models we 
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expect Uc ~ -0/(1 + J/U), decreasing with J/U. But for 
J » A, the configuration is the 5 = 2 high-spin state, 
and depends on the electron filhng of the high-spin 
state. In the four-orbital model, the system is at half- 
filling, ^ U + 3J, and Uc - D/{1 + 3J/U), still 
decreasing with J/U. But in the five-orbital model with 
n = Q, G^ ^ U - 3J. This gives Uc ~ D/{1 - 3J/U), 
which increases with J/U. Therefore, we expect a non- 
monotonic behavior of Uc with increasing J/U, which is 
seen in the numerical results shown in Fig. [1] 

The Hund's rule coupling also affects the properties in 
the metallic state. In Fig. [5J we compare the evolution 
of quasiparticle spectral weight Za for the same model 
but at two different J/U values. In both cases, the insu- 
lating phase is the 5 = 2 high-spin Mott state. On the 
other hand, Za behaves very differently in the metallic 
phases. At J/U = 0.1, Uc ^ D, and Za drops rapidly 
down to zero only at U ~ Uc- The orbital dependence 
of Za is weak. At J/U = 0.25, the Mott transition takes 
place at U ^ 3D. But Z^ drops rapidly to small but 
non-zero values at f7 w 2.7 cV. This rapid drop allows us 
to defines a crossover scale U* in the metallic state. In 
Fig. [TJ we plot this crossover line in the phase diagram. 
For large J/U, U* can be smaller than D. U* increases 
with decreasing J/U and the crossover line ends when 
it crosses the MIT phase boundary at J/U « 0.11. At 
U < U* the spectral properties of the system is similar to 
its non-interacting limit, with weakly renormalized quasi- 
particle spectral weights. But for U > U* , the quasipar- 
ticle spectral weights are strongly suppressed. In this 
regime, electron correlations are sufficiently strong in the 
metallic phase even at C/ < £>, and the system is close to 
a Mott insulator. So U* roughly separates the regimes 
of a weakly correlated metal and a strongly correlated 
metal. The strongly correlated metallic phase has the 
features prescribed in the incipient Mott picture. As an- 
other remarkable observation, we find that for U > U* , 
Za become significantly orbital dependent. For exam- 
ple, at J/U = 0.25 and for [/ > 4 cV, Z^^y < 0.1 and is 
much smaller than that of the other orbitals. This im- 
plies that the system is close to an orbital selective Mott 
phase (OSMP). We have calculated the groundstate en- 
ergy of an OSMP with xy orbital insulating but all other 
orbitals metallic, and compare it with the energies of the 



metallic and insulating solutions in the inset of Fig. [2] 
Though the OSMP never becomes the true ground state, 
it is indeed energetically close. The stabilization of an 
OSMP requires a high-spin configuration^. This is con- 
sistent with the observation that a threshold is needed to 
trigger the strongly correlated metallic state. 



IV. CONCLUSION 

We have developed a [/(I) slave-spin theory, which al- 
lows the study of metal-to-Mott-insulator transition in 
both single- and multi-orbital systems. For models with 
a single orbital or multiple degenerate orbitals, we show 
that the mean-field theory in the slave-spin-condcnsed 
phase is mathematically equivalent to that of the pre- 
vious Z2 slave-spin mean-field theory. For models with 
multiple non-degenerate orbitals, our C/(l) slave-spin for- 
mulation provides proper descriptions for both the metal 
and Mott-insulating phases. 

We have applied the U{1) slave-spin approach to study 
a five-orbital Hubbard model for the parent iron pnic- 
tides. We find that the model exhibits a metal-to-Mott- 
insulator transition. The interplay between the Hund's 
rule coupling and crystal field splittings strongly affect 
the Mott transition and the associated phases. The insu- 
lating phase can be either an 5* = 1 low-spin Mott state or 
an S* = 2 high-spin Mott state, depending on the strength 
of the Hund's rule coupling. In the metallic phase, a 
crossover between a weakly correlated to a strongly cor- 
related metallic phase exists when the Hund's rule cou- 
pling is beyond a threshold. Inside the strongly corre- 
lated metallic phase, the quasiparticle spectral weights 
are strongly suppressed, in agreement with the incipi- 
ent Mott picture. In this phase and in the vicinity of 
the Mott transition, we find that an orbital selective 
Mott phase has ground state energies which are nearly as 
competitive as those of the metallic and Mott insulating 
states. 
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